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ABSTRACT 

Deep F555W and F814W Hubble Space Telescope ACS images are the basis for 
a study of the present day mass function (PDMF) of NGC 346, the largest active 
star forming region in the Small Magellanic Cloud (SMC). We find a PDMF slope 
of T = —1.43 ± 0.18 in the mass range 0.8 — 60 M , in excellent agreement with 
the Salpeter Initial Mass Function (IMF) in the solar neighborhood. Caveats 
on the conversion of the PDMF to the IMF are discussed. The PDMF slope 
chang GS, clS 8b function of the radial distance from the center of the NGC 346 
star cluster, indicating a segregation of the most massive stars. This segregation 
is likely primordial considering the young age (~ 3 Myr) of NGC 346, and its 
clumpy structure which suggests that the cluster has likely not had sufficient 
time to relax. Comparing our results for NGC 346 with those derived for other 
star clusters in the SMC and the Milky Way (MW), we conclude that, while the 
star formation process might depend on the local cloud conditions, the IMF does 
not seem to be affected by general environmental effects such as galaxy type, 
metallicity, and dust content. 
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INTRODUCTION 



Since the pioneering study by iSalpeterl (119551 ). one of the major issues in astrophysics 
is whether the initial mass function (IMF), the relationship that specifies the mass distri- 
bution of a newly formed stellar population, is universal or, alternatively, determined by 
environmental effects. In particular, establishing whether the IMF has been constant over 
the evolution of the universe, or if it varies with time (redshift) and/or metallicity has crucial 
consequences on the evolution, surface brightness, chemical enrichment and baryonic content 
of galaxies, and on the evolution of light and matter in the universe. 

Due to a number of effects and observat ional uncertainti es such as mass segregation, field 
star contamination, and variable extinction (jElmegreen!l2006l ) . the "true" IMF is very difficult 
to measure. Despite these limitations, star clusters still provide the best determination 
of the IMF, because they are essentially single-age, single-metallicity stellar systems. In 
gener al, all star clusters lose more than 80% of their stellar content in their first ~ 10 
Myr (IKroupa &: Boilyll2002l ; lLada fc Ladall2003l ). Therefore, in order to perform an accurate 
census of the initial stellar content it is necessary to focus on the youngest stellar systems. In 
addition, young (age < 3 — 5 Myr) star clusters in nearby galaxies (e.g. the Large -LMC- and 
Small Magellanic Cloud -SMC) offer to us the unique advantage to individually resolve - 
because of their proximity — statistically significant samples of stars down to the sub-solar 
mass regime. 

In this paper, we foc us our attention on the mass functi on (MF) of NGC 346, an 
extremely young (~ 3 Myr, iBouret et al.ll2003l : ISabbi et al.l 120071 . hereafter S07), and active 
star cluster that excites the largest and brightest Hn region (N66) in the SMC. NGC 346 



K-2000 = 00 fe 59 m 5.2-,^,nnri 



in the entire SMC ( Walborn 



— 72 deg 10'28") contains a maior fraction of the O stars known 



19781 : IWalborn fc Bladeslll986l : iNiemela et alll986l : iMassev et al. 



19891 ). The bright end of its stellar population (V < 19.5 mag) has been well investigated in 
the past 20 years. Spectral inve stigations of the brightest members identified several stars 



of spectral type 6 .5 or earlier ( Massey et al.l Il989l ; iHeap et al.l 120061 ; iMokiem et al.l 12006 



(|l989|) found the MF for the brightest stars in NGC 346 



Evans et al.ll2006l ). IMassev et al. 
- down to 5 M — to have a slope that is consistent with what has been found for massive 
stars near the Sun and in the LMC. 

High resolution observations obtained with the Advanced Camera for Survey (ACS) 
on board of the Hubble Space Telescope (HST) recently revealed that the star formation 



1 Based on observations with the NASA/ESA Hubble Space Telescope, obtained at the Space Telescope 
Science Institute, which is operated by AURA Inc., under NASA contract NAS 5-26555 
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history of the region where NGC 346 is located is quite complex. There is evidence for 
old episodes of star formation (SF), in the field of the SMC approximately between 3 and 
5 Gyr ago. After that, the SF activity in this region appears to have been significantly 
diminished, with a possible moderate enhancement ~ 150 Myr ago (S07). An intermediate- 
age (~ 4.5 Gyr) stellar cluster (BS90), characterized by a core radius r c ~ 25arcsec and a 
tidal radius r t ~ 130 arcsec, is located at a projected distance of ~ 23 pc northeast from 
the center of NGC 346 (S07). Spectral analysis of NGC 346 members revealed the presence 
of stars as young as ~ 1 Myr, (IMassey et al.ll2005l ). but stars ~ 5 Myr old are also present 
Jrleap et al.lbood : iMokiem et alil2006h . 



Star formation is likely still active in the region: hu ndreds of pre-mai n sequence (pre- 
MS) stars in the mass range 0.6 — 3 M were discovered by lNota et al.l (120061 ) from HST / ACS 
images. Spitzer Space Telescope (SST) observations also detecte d a myriad of emb edded 
young stellar objects (YS Os) scattere d acro ss the entire region (IBolatto et al.l 120071 ) . By 
assuming a Salpeter IMF, ISimon et al.l (120071 ) calculated that more than 3000 M Q have been 
formed in the last ~ 10 6 yr, concluding that > 6% of the current SF in the SMC is taking 
place in NGC 346. 

This recent stellar population does not appear to be uniformly distributed within the 
ionized nebula, but is rather organized in many, likely coev al, subclusters (S07) t hat are coin- 



cident with clumps of molecular gas previously detected by lRubio et al.l (120001 ) . ISimon et al. 
(120071 ) noted that almost all the subclusters host one or more YSOs. 



The high sensitivity of the HST /ACS allowed us to investigate the stellar content of 
NGC 346 from ~ 6OM down to ~ O.6M , making NGC 346 one of the few known regions 
where a MF can be determined over two orders of magnitude (th e other classical case being 
the Orion trapezium cluster in our own MW; lElmegreenl 120061 ) . This paper gives a new 
derivation of the NGC 346 MF, reviews the strengths and limitations of the methods, and 
presents the conclusions that we can conservatively derive. The paper is organized as follows: 
a short description of the data reduction procedure is presented in $2) in $3] we present the 
color-magnitude diagram (CMD), and we describe the stellar populations identified. In Sgl 
we present the MF we obtained for the NGC 346 cluster, we discuss its spatial variations, 
and we analyze the impact of the environment on the MF. The results of this paper are 
discussed in 351 
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OBSERVATIONS AND DATA REDUCTION 



Multiple, deep, images of NGC 346 were obtained during Cycle 14 using the Wide 
Field Channel (WFC) of the HST/ ACS (GO-10248; PI A. Nota). A detailed description of 
the observations and of the stellar photometry has been presented in S07. Several exposures 
were taken through the filters F555W (~V) and F814W (~I). A dither pattern was especially 
designed to allow for hot pixel removal, improve the sampling of the point spread function 
and the accuracy of the photometry by averaging the flat-field errors and by smoothing 
over the spatial variations in the detector response. Short and long exposures were taken to 
maximize the photometric range. 

The entire dataset was processed adopting the standard STScI ACS calibration pipeline 
(CALACS), and exposu res in each band were co-added using the MULTIDRIZZLE package 
( IKoekemoer et al.ll2002l ) to correct for the geometrical distortions, and remove cosmic rays 
and hot pixels. The final exposure time of the deep combined F555W and F814W images is 
~ 4100 s, and the total area imaged i s ~ 5' x 5', correspo nding to ~ 88 x 88 pc, assuming 
a distance to the SMC of 60.6 Kpc (IHilditch et al.ll2005l ). The photometric reduction of 
the images was performed with the DAOPHOT pack age within the IRAF §| environment. 
We followed the sam e methodology applied to R136 ( Sirianni et al. 200ol ) and NGC 330 
( iSirianni et al.ll2002l ) and used the PSF fitting and aperture photometry routines provided 
within DAOPHOT to derive accurate photometry of all stars in the field (see S07 for more 
details). 

The instrumental magnitudes were corrected for the charg e transfer efficiency ( CTE) 
effects that affect ACS data, using the parametric correction from lRiess fc Mack! (120041 ) . The 
photometric calibration to the Vegamag system was performed by converting the magnitudes 
of the individ u al sta rs to an aperture radius of 0.5" and applying the zero points listed in 
Sirianni et al.l (120051 ). 



A SMC field (a = 00:58:42.5; 5 = -72:19:46) was observed with WFC for comparison. 
The observed field is located ~ 9'16" south (~ 163 pc in projection) from NGC 346 (Fig. [T]). 
The data were analyzed following the same procedure described above. 

Extensive artificial star tests were designed to accurately assess the level of completeness 
of the photometric data. More than 3,000,000 artificial stars — obtained from the scaled 
PSF used for the photometric analysis of the frames — were added to both the F555W and 
the F814W images using the subroutine addstar in DAOPHOT, and subsequently retrieved. 



2 IRAF is distributed by the National Optical Astronomy Observatory, which is operated by AURA Inc., 
under cooperative agreement with the National Science Foundation. 
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The artificial stars were distributed onto the F555W image randomly in magnitude using 
a step function. This function extends two magnitudes below the detection limit of our 
observations in order to probe with sufficient statistics the range of magnitude in which 
incompleteness is expected to be most severe. This is a crucial step to assess the robustness 
of the MF. 

In each run ~ 7000 artificial stars were distributed on a regular grid to avoid an unreal- 
istic increase of the crowding conditions in the analyzed frame. Then, the frame to which the 
artificial stars were added was processed exactly as the original frame. We determined the 
level of completeness of the photometry by comparing the list of artificial stars added with 
the list of the stars recovered. In our experiments an artificial star is considered recovered if 
it is found in both F555W and F814W frames with an input-output difference in magnitude 
Am < 0.75, and at the same time satisfying all the photometric selection criteria. 

The artificial star experiments show that the completeness is a function of both magni- 
tude and color: our photometry is complete at a 50% level down to magnitude mF8i4W — 25 
at color m F 555w — m F8uw — —0.2. At redder colors (e.g. m F 555w — m F8i4W — 2.0) the 
photometry is complete at the 50% level only at magnitude m F si4w — 23.2 (Fig. [2]). The 
50% level of completeness defines the magnitude limit we have conservatively imposed for 
the MF determination. 



3. THE STELLAR POPULATIONS IN NGC 346 

Figure [2] shows the mpsi4w vs - m F555W — m F8uw CMD of all the stars identified in the 
WFC NGC 346 data with (Jdao < 0.1 mag. We can distinguish in this CMD the different 
stellar populations that coexist in the region. 

• A rich old stellar population, with a main sequence (MS) extending from mpsuw — 
22.0 down to mpsuw — 25.5. The evolved stars of this population clearly define the red giant 
branch (RGB), with the brightest stars at m F814W ~ 15.0, and a tight red clump (RC) at 
^F8i4W — 18.5. The narrow subgiant branch (SGB), visible at magnitude mp8i4w — 21.0, in 
the color range 0.45 < mpsssw — ^F8i4W < 0.95, indicates that, in the NGC 346 region, the 
old stellar population is dominated by the presence of the BS90 star cluster. The presence 
of SGB stars both fainter and brighter than the BS90 SGB indicates that, in the field, SF 
occurred before and after the formation of this star cluster. 

• A young population. The bright (12.5 < mp8i4w < 22.0) and blue (—0.3 < mpsssw ~ 
^fsuw < 0.4) MS, well visible to the upper left of the CMD, is the most remarkable feature 
of the youngest stellar population. The majority of these stars belong to the cluster NGC 346, 
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and have an age of ~ 3 Myr. 

• Hundreds of pre-MS stars are visible between 0.6 < m F5 5 5W — m F 8i4w < 2.2, below 
magnitude m F814W ~ 22.0. 



THE MASS FUNCTION 



The MF, commonly indicated by £, is defined as the fractional number of stars per mass 
interv al per unit area. Among the various parameterizations of the MF (see iKroupa &: Boily 
20021 ) . one of the most used is the one proposed by IScalol (119861 ). where the MF is char- 
acterized by the logarithmic derivative T = rf(log^(log(m))/ dlogjTt 3 where £ (log(m) is the 
MF, and T is its slope. The reference IMF, as derived by ISalpeterl (119551 ) for the solar 
neighborhood, has a slope T = —1.35. 

The MF of a system is obtained by counting the number of stars in mass intervals. 
Under the assumption that all the stars of a system formed at the same time, the MF can be 
derived by counting the number of stars as a function of magnitude (luminosity function), 
which translates into the n umber of stars per un it stellar mass by using an appropriate mass- 



luminosity relation (see e.g. lde Grijs et al. 



2002J). However, this approach can be applied only 



to those stars currently on the MS, where the relationship between the luminosity of a star 
and its mass is reasonably well known, and assumes that all the stars have the same age, 
and the same metallicity, or, in other words, that we are studying a simple stellar population 
(SSP). 

As discussed earlier, there is evidence that NGC 346 formed stars over the last 5 Myr. 
Furthermore hundreds of stars in NGC 346 are still in the pre-MS ph ase (where their lu- 
minosity strongly depends on age), and star formation is still ongoing (IBolatto et al.l 120071 ; 
Simon et al.l 120071 ). For all these reasons we cannot model NGC 346 as a perfect SSP. To 
overcome the age-luminosity degeneracy wh ich comes from the fact that NGC 346 is not a 
SSP, we applied the technique presented by iTarrabl (119821 ) and iMassey et al.l (119951 ). which 
consists of directly counting the stars between two theoretical evolutionary tracks, accord- 
ing to their position in the Herztprung-Russell Diagram (HRD). This technique has the 
advantage of avoiding any assumptions on the age of the stellar population. 

Because in NGC 346 SF occu rred over an interval of time comparable with its age and 
is still ongoing (jSimon et al.ll2007l ). it is almost impossible to derive its IMF without making 
assumptions on how the star formation rate varies in time, and on the mass distribution of 
the stars which are currently forming. For this reason we will focus our attention on the 
present-day mass function (PDMF) of those stars that are already visible in the optical 
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bands, but have not yet evolved off the MS. 

To infer the PDMF of NGC 346, we used two different sets of e volutionary tracks : 



for stars below 7 M , we adopted pre-MS stars evolutionary t racks by ISiess et al.l (120001 ) 



while we used Padua evolutionary tracks (IFagotto et al.l Il994l ) for stars above 9 M . In 



order to minimize the uncertainties due to the photometric conversions, we transform the 
evolutionary tracks into the Vegamag system by applying the convertions calculated by 



Origlia fc Leithererl (120001 ) . The adopted evolutionary tracks, superimposed on the m F555W 



vs. m F555W — m F8UW CMD, are shown in Figure [2J 

We computed the PDMF between 0.8 and 60 M using the following binning: 0.8, 1.2, 
1.6, 2.0, 2.5, 3.0, 4.0, 5.0, 6.0, 7.0, 9.0, 12.0, 15.0, 20.0, 30.0, and 60.0 M . 

To properly account for the effects of crowding and incompleteness, which varies with 
both magnitude and color (e.g. the dashed line in Fig. [2] shows the 50% level of completeness 
of our photometry), stars identified between two evolutionary tracks were divided into bins 
of 0.5 in magnitude and 0.25 in color. The number of stars identified in each bin was then 
normalized to the logarithmic width of the mass range spanned by the tracks and to the area 
of the observed region. We then applied the appropriate completeness factor. Figure [3]-panel 
(a) shows the PDMF obtained with this approach. 

The resulting PDMF reflects the composite stellar population in this region. In order to 
derive the "true" PDMF for the NGC 346 cluster, we need to account for the contamination 
of the old cluster BS90 and of the SMC field. 

• The BS90 old cluster. We used the King model that best fits the stellar density profile 
of BS90 (S07) to determine, at each spatial position, the expected stellar contribution from 
BS90. We selected the stars that lie within 20 arcsec from the projected center of gravity 
of BS90 as representative of the old cluster stellar population. We then scaled the number 
of stars to the estimated stellar density, and removed them from our PDMF. We found that 
~11% of the stars are likely to belong to BS90. Figure [3]-panel (b) shows the PDMF after 
the subtraction of the contribution from BS90. 

• The SMC field. We used a SMC field (a J20 oo = 00 /l 58 m 42.5 s , 5 J20 oo = -72°19'46"), 
located at a projected distance of ~ 163 pc from NGC 346, observed with WFC in the 
same filters, to remove the contribution of the SMC stellar population from our data. This 
contribution is important, given that the surface stellar density in this region is ~ 9.8 stars 
pc~ 2 . As we did for NGC 346, we counted the number of stars lying within two evolutionary 
tracks, normalized to the logarithmic width of the mass range spanned by the tracks and to 
the area of the observed region, and then we applied the appropriate completeness factor. 
Then for each mass bin, we subtracted the estimated number of stars belonging to the field 
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from the PDMF of NGC 346. From the comparison between the reference field and NGC 346 
region, we can conclude that in NGC 346 region about 54% of the stars below 12 M belong 
to the SMC field. 

The final PDMF of NGC 346, after the subtraction of the SMC field and BS90, is shown 
in Figure [3]-panel (c). A weighted least mean square fits of the data indicates that, between 
0.8 and 60 M , the slope of NGC 346 PDMF is T = -1.43 ± 0.18, in good agreement with 
the value derived by Salpet er for the IMF of th e solar neighborhood. Furthermore we note 
that, as already derived by iMassey et al.l (119891 ). the PDMF is quite steep above 5 M (we 
derive a slope T = 1.87 ± 0.41), and becomes flatter below this value of mass. 

Some uncertainties still affect the derived slope of the PDMF: 



• As IMassey et al.l (119951 ) have shown, optical and near-UV colors only sample the low- 
frequency tail of the spectral energy distribution of the most massive and hot stars; as 
a consequence, optical photometry is not adequate to estimate the PDMF of massive 
stars, and spectroscopy is nec essary. In Figure EH c) we have compared our derived 
PDMF with that published by IMassey et al.l (119891 ) (open triangles): we found a very 
good agreement in the overlapping range of masses below 20 M . As expected, on the 
basis of photometry alone, we underestimate the masses of the most massive stars. 

• The presented PDMF does not take into account uncertainties due to unresolved binary 
systems, which may be important. In the O rion Nebula Cluster, f or example, each 
massiv e star has, on average, 1.5 companions (IPriebisch et al.lll999l ). ISagar fc Richtler 
(Il99ll ) estimated that, if each star in the mass range 2 to 14 M has one companion, 
the average IMF slope they derived for five young clusters in the LMC would change 
from T = —1.3 to —1.7. This correction, however, depends on t he distribution of th e 
companions, which is not yet well-known (IPriebisch et al.lll999l : iDuchene et al.ll200ll ). 
A spectroscopic survey of NGC 346, carried on with the Fiber Large Array Multi- 
Element Spectrograph (FLAMES) instrument at the Very Large Telescope (VLT), 
indicates a binary fraction of ~ 26% for the most massive stars (5 < M < 60 M , 



Evans et al.ll2006l ). 



• The derivation of the PDMF requires an accurate knowledge of the amount of extinc- 
tion. The foreground extinction to the SMC is low, and the only appreciable presence 
of gas and dust along the line of sight of NGC 346 is likely associated with the star- 
forming region itself. A visual inspection of the region shows that dust and gas are not 
uniformly distributed, but detailed reddening maps for the region are not yet available. 



• The real slope of the MF can be altered by an inaccurate subtraction of contaminating 
field stars: dust and gas for example can obscure field stars, causing an over subtraction 
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of the low-mass stars, and a flatter MF. Local differences in the stellar density or in 
the distance dispersion between the analyzed region and the reference field may also 
artificially alter the shape of the MF (see §4. ID . 



A major source of uncertainty comes from the pre-MS evolutionary t racks, which do 



not yet well represent o bservations (see for example the discussions in IChabrierl 12003 



Hillenbrand et all 120071 ). 



4.1. The spatial variations of the PDMF 

The high spatial resolution and the depth of our photometric data have allowed us to 
study in detail the spatial variations of the PDMF within NGC 346. Already a cursory 
inspection of the image shows that the brightest stars are found in the central region. Does 
this imply that the most massive stars are concentrated in the cluster center? 

To answer this question, we have calculated the variation of the mass function with 
projected distance from the NGC 346 center. In order to determine the position of the 
projected center of gravity of NGC 346 we first selected the stars belonging to the young 
MS (defined by magnitude m F si4w brighter than 21.0 and colors m F555W — m F814W bluer 
than 0.4). Then we drew a circle with a radius of 50 pixel around each star and counted 
the number of stars falling within each circle. We computed the weighted average of the 
MS stars a and 5 coordinates, using the sum of the stars within each circle as weight. The 
position we obtained for the stellar association is a = 00 /l 59 m 05'/8 5 = -72°10 m 35", [J2000], 
with a la uncertainty in both a and in 5 of ~ 1, that corresponds to about 25 pixels in the 
HST/WFC images. 

To study the spatial variations of the PDMF, we divided the region in four annuli, at 
radial distances from the center R< 4.00, 9.00 (which corresponds to the half mass radius), 
14.00, and 19.80 pc (Fig. 0J. Less than 2.6 stars/pc 2 belong to NGC 346 outside ~ 20 pc 
from the center. PDMFs derived for annuli 1, 2, 3, and 4 are shown in Figure and the 
corresponding slopes, as a function of the distance from the center, are shown in Fig [6j 
In the innermost region the PDMF appears quite flat, with T = —1.03 ± 0.14 (Fig 
It becomes steeper moving from the center to the periphery, where its exponent becomes 
r = -2.08 ±0.14. 

We applied a Kolmogorov-Smirnov (K— S) test to verify the significance of the differences 
found between the PDMF obtained in the innermost annulus (annulus # 1) and those derived 
for the other three annuli. According to the K-S test, the probability that the PDMFs 
derived in annuli 2, 3, and 4 are drawn from the parent distribution of the PDMF in the 
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innermost annulus is 3.9%, 0.3%, and 0.06% respectively. Therefore, we can safely exclude 
that regions 3 and 4 are drawn from the same parent distribution of region 1. 

The projected density of massive stars from the innermost annulus to the periphery 
decreases by a factor of 60, whereas low-mass stars are depleted only by a factor of 6, 
indicating that the change in the PDMF slope is due to a lack of massive stars in the 
periphery, rather than an excess of low mass stars there. A similar trend was observed 
in Or ion Nebula Cluster (ONC). Discussing the excess of massive stars in the center of the 
ONC Hillenbrand! (119971 ) pointed out that mass distribution could be biased by an age effect, 
with the most massive stars also be ing the younge s t. Th is does not seem to be the case for 
NGC 346: spectroscopic analysis by lMokiem et al. J2006h have demonstrated that, even if a 
spread of about 5 Myr in the ages of NGC 346 stars is present, stars of different ages appear 
uniformly distributed over the whole cluster, and SST infrared observations indicate that 
there is ongoing SF throughout the com plex with a lower limit on the star formation rate of 
3.2 x lO^MQyr- 1 Jsimon et alihoOTT ). 



As already mentioned in §U dust and gas can obscure field stars, and this may cause an 
over subtraction of the low-mass stars, resulting in a flatter MF. To exclude the possibility 
that the spatial variation of the MF can be totally ascribed to variations in the gas and 
dust densities, we derived the MF of the innermost region also under the assumption that 
the contamination from the SMC field is negligible, and in this case we obtained a value 
for T = —1.11 ±0.13. We also estimated the surface stellar density necessary to obtain 
the same MF slope in the innermost region that we find for the outermost one. In order 
to obtain a slope of T = —2.08 we should observe a stellar density of ~ 1000 star pc~ 2 
between 0.8 — 1.2 M , which is almost a factor of 2 higher than the observed value; such a 
high stellar density makes it unlikely that the derived variation can be ascribed only to a 
poor subtraction of the SMC field. 



4.2. Is the mass segregation primordial? 

As for many other young star clusters, NGC 346 shows a shallower PDMF in its central 
region. It is still a matter of debate whether the segregation of the most massive stars in 
young clusters is due to evolutionary effects, or instead is due to the initial conditions of the 
cluster formation. In the first case the most massive stars form elsewhere in the cluster and 
then sink to the center through dynamical interactions with the numerous low-mass stars. 
On the contrary, if the mass segregation is due to initial conditions, then it reflects how the 
cluster formed. 
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As discussed by S07, NGC 346 shows a very complex morphology, with the stars as- 
sembled in many coeval sub-clusters. Due to the lack of information on the dynamics of the 
sub-clusters, at zero order we can consider all the stars in NGC 346 as members of an indi- 
vidual star cluster. Thus, to understand whether the mass segregation found in NGC 346 is 
primordial or not, we can compare its age with its mass segregation timescale (T msg ), defined 
as T msg = 2T rc i ax m av /m max , T re i ax is the relaxation time of t he cluster, m a v is the average 



stellar mass and m max is the mass of the most massive star (jKroupal 120041 ) . For NGC 346 



we found T msg ~ 2 x 10 7 yr (values used to derive the T msg are listed in Table [Q, one or- 
der of magnitude larger than the age of NGC 346, supporting the idea that the observed 
segregation of the massive stars is likely due to initial conditions, rather than dynamical 
evolution. 

It has been suggested that NGC 346 is probably the result of the collapse and subse- 
quent fragmentation of the initial giant molecular cloud into multiple "seeds" of SF. The 
fact that the observed subclusters appear almost coeval, and that many of them appear 
to be connected by arcs of dust and gas, seems to be a good observational match to the 
conditions predicted by the hierarchical fragmentation of a turbulen t molecular cloud model 
(klessen fc Burked 1200(1 iBonnell k Batelbooi iBonnell et al1l2003h . If this is the case, the 
fact that the various "seeds" of star formation are still distinguishable would imply that 
NGC 346 is not completely relaxed, and that the majority of the stars would be still close 
to their birth positions. Also in this case the observed mass segregation would be primordial 
and would represent a feature of the way the cluster is forming. 



4.3. Environment and the IMF 



The Magellanic Clouds differ from the Milky Way in a variety of ways, including their 
mass densities, tidal fields, chemical abundances, and experiences with recent external per- 
turbations. Both the Magellanic Clouds also are well known for their large populations of 
relatively massive star clusters: densities of host galaxies and the densities of their star form- 
ing sites evidently are not simply related. In the SMC the star formation process on large 
spatial scales itself differs from that in the Galaxy, more li kely being a result of cloud-cloud 
interactions (jZaritsky et al.ll2000l; IZaritsky fe Harris! 120041 ) and the turbulent generation of 



self-gravitating structures (e.g. lElmegreenl 120001 ) . In some regions this may be further en- 
hanced thr ough the formation of dens e shells and supershells which are sites for further star 
formation (IHatzidimitriou et al.l 120051 ). Thus the SMC, with its low metallicity and exter- 
nally driven star formation in a gas-rich system, provides insights into the development of 
stellar populations in small pre-galactic systems. 
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The form of the stellar IMF in the SMC has implications for a wide range of studies. 
The IMF has been studied for higher mass stars in the SMC: the derived slopes are in 
good agreement wi t h the values ob tained, in similar ranges of masses, in the MW (e.g. 
Massey et al.l Il989l ; iHill et al.l 119941 ) . Simila rly, existing investigations of the IM F down 
to 1 M yield Salpeter-like IMF slopes (e.g. iMateol Il988l : Ichiosi fc VallenarilboOTT ). Thus 
while the star formation process might be somehow affected by the conditions of the local 
interstellar medium, evidently the change in conditions between the SMC and spirals has 
not significantly affected the IMF. 

Most of the stellar mass resides in stars with M<1 M . Thus our extension of the lower 
limit for the PDMF in NGC 346 to 0.8 M Q in this study adds further support to the overall 
normalcy of the IMF. Since this result comes from one major star forming complex, our data 
further suggest that this "normal" IMF is imprinted during major formation processes. 

On a smaller spatial scale, environm ent, however, d oes a ppear to matter within the 
NGC 346 complex. Consistent with the ISirianni et al.l (120021 ) work on NGC 330 in the 
SMC we find that while the overall PDMF in NGC 346 follows a Salpeter form, there are 
significant local variations. The NGC 346 complex offers, therefore, an ideal opportunity to 
characterize the occurrence and the nature of the local variations. This will be the objective 
of future work. 



5. DISCUSSION AND CONCLUSIONS 



Our analysis of the stellar content of NGC 346 indicates that, within all the uncertainties 
that can affect the determination of a MF, the PDMF slope o f NGC 346, at least down to 
the solar mass regime, is consistent with the value derived by ISalpeterl (119551 ) for the IMF 
in the solar neighborhood, further confirming that the stellar MF does not strongly depend 
on the metallicity, or the mass or morphological type of the parent galaxy. 

We studied how the PDMF varies as a function of the distance from the center of 
NGC 346. We found that the PDMF is quite flat in th e center, and becomes steeper in the pe- 



riphery. As already found in other young star clusters (iHiilenbrandl 199 71 ; Hillenbrand fc Hartmann 



19981 ; ISirianni et al.l 120021 ) the increase in the MF slope with the distance from the center is 
due to the lack of massive stars in the periphery rather than an excess of low-mass objects 
there. The projected density of massive stars from the innermost annulus and the outskirts 
decreases by a factor of 60, whereas low-mass stars are depleted only by a factor of 6; the 
most massive stars are segregated in the center. 



We conclude that the spatial stellar mass distribution is unlikely to be biased by age 
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effects. SST observations indicate that there is ongoing SF throughou t the complex wit h 



a lower limit on the total s tar fo rmation rate of 3.2 x 10 3 M Q yr 1 (Isimon et al. 2007 ). 



Furthermore iMokiem et al.l (120061 ) noted a spread of about 5 Myr in the ages of NGC 346 
stars, but did not find any correlation between the age and the position of the stars over the 
whole cluster. 

S07 noted that stars in NGC 346 appear to be organized in several coeval sub-c l usters , 
embedded in Hll gas, and coinciding with the CO clumps analyzed bv lRubio et al.l (120001). 
YSOs between Class I and Class III are found in 14 of the 16 identified sub-cluster ( Simon et al 



20071 ). The sub-clusters are also connected by filaments and arcs of gas and dust. The com- 
plex structure of NGC 346 and its young age, compared to its dynamical time, suggest that 
NGC 346 is not dynamically evolved, and probably not yet relaxed. On the basis of these 
considerations we conclude that the observed segregation of the most massive stars is likely 
primordial, and reflects how the cluster formed. 

S07 suggested that in NGC 346 most of the star formation was likely not triggered, but 
rather resulted from the turbulence driven density variations within a giant interstellar cloud 
complex, con ditions predicted by the hierarchical fragmentation of a turbu l ent molecular 



cloud model (lElmegreenll2000l ; iKlessen fc Burkertll2000l ; Bonnell fc Batdl2002l : iBonnell et al. 



20031 ). According to this model, the fragmentation of the cloud is due to supersonic turbulent 
motions present in the gas. The turbulence induces the formatio n of shocks in the gas, and 
produces filamentary structures (IBate. Bonnell. &: Brommll2003l ). The chaotic nature of the 
turbulence increases locally the density in the filamentary structures. When regions of high 
density become self-gravitating, they start to collapse to form stars. Simulati ons show that 



star f ormation occurs simultaneously at several different locations in the cloud (IBonnell et al. 
20031 ). as appears to be the case for NGC 346. 



In this scenario, stars forming near t he cluster center would haye hi gher accretion 
rat es, due to th e loca l higher gas density (ILarsonl Il982l ; IBonnell et alll200ll ). Simulations 
by IBonnell et all (120011 ) indicate that although low-mass stars form equally well throughout 
the entire cluster, the most massive stars are almost exclusively segregated to the central 
regions. It is interesting to note that the two most massive stars (one of them being an 03V 
spectral type, Figs. |7|) are ~ 7pc off the center of the cluster (outside the core of NGC 346, 
but still within the half mass radius). Did these stars form where they are observed today, 
in a quite low stellar density region? 



Vine fc Bonnelll (120031 ) studied the effect of strong stellar winds and photoionizing fluxes 



coming from O- and B-type stars on the gas content in a young star cluster. According to 
their simulations, the gas expulsion is unlikely to remove any initial mass segregation, but 
it can still result in the ejection of some of the most massive stars to positions outside the 
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cluster core, although still within the half mass radius. According to iBouret et al.l (120031 ) 
mechanical power from O-star winds is not a dominant factor in the evolution within N66, but 
the molecular g as in the original c loud has been strongly photodissociated and photoionized 
by the O stars Jllubio et aliboool ). 



Niemela et al.l (119861 ) measured the radial velocity of five of the innermost and brightest 
stars in NGC 346, and found a velocity dispersion of ~ 5kms _1 , indicating that, even in 
~ 1 Myr, they can be formed in the center of the cluster and then ejected where they are 
observed today. Another possibility is that these stars formed where are ob served, as the 
resul t of a very efficient competitive accretion process into a potential well (IBonnell et al. 
2003h . 



Further spectroscopic analysis of the dynamics of the gas and the sub-clusters can 
confirm whether the giant molecular cloud which formed NGC 346 underwent under a hier- 
archical fragmentation. The initial superson ic turbulence originates shocks in the gas which 
rapidly remove kinetic energy from the gas (lOstriker. Stone, fc Gammidl200ll ): in this case 
we expect to find low gas velocities around the sub-clusters. A detailed study of the sub- 
clusters dynamics will also reveal if NGC 346 is or is not already relaxed and virialized. We 
will present these results in a future paper. 
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Table 1: NGC 346 properties derived in this paper. 



NGC 346 MF slope between 1 and 60 M© T = -1.43 ± 0.18 

MF slope within 4 pc from the center T = —1.03 ± 0.14 

MF slope between 4 and 9 pc from the center T = —1.36 ± 0.09 

MF slope between 9 and 14 pc from the center T = -1.43 ± 0.09 

MF slope between 14 and 20 pc from the center T = -2.08 ± 0.14 

NGC 346 half mass radius 9 pc 

NGC 346 total mass M Tot ^ 3.9 x 10 5 M 

NGC 346 crossing time T cross ~ 1.8 x 10 6 yr 

NGC 346 relaxation time T rclax > 5.7 x 10 8 yr 

NGC 346 mass segregation time T msg > 2 x 10 7 yr 
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Fig. 1. — Relative position of NGC 346 and of the field observed as reference (a). An 
enlarged image of NGC 346 is shown in panel (b), together with an enlarged view of the 
field (c). 
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Fig. 2. — NGC 346 m F814H / vs. m F555W — m F814W / CMD. We show Padua MS evolutionary 
tracks for stars more massive than 9 M (grey dotted lines) and Siess Pre-MS evolutionary 
track for stars in the mass range between 0.8 and 7 M (grey continuous lines). The grey 
dashed line indicates the 50% level of completeness of our photometry, as derived from the 
artificial star experiments. 
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Fig. 3. — PDMF of NGC 346 in the mass range between 0.8 and 60 M . Panel (a) shows 
the derived PDMF before removing the contributions from the old star cluster BS90 and the 
field of the SMC; panel (b) shows the PDMF after the correction for the cluster BS90, and 
panel (c) the "true" PDMF of NGC 346, after the su btraction of BS90 an d the field of the 
SMC. Open triangles represent the IMF as derived by iMassey et al.l (119891 ) . 




Fig. 4. — NGC 346 WFC/ACS mosaiked data. We show the four annuli used to study the 
spatial variations of the PDMF. 
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Fig. 5. — PDMFs of the four regions of NGC 346 analyzed: in panel (1) we present the 
PDMF derived within 4 pc from the center of the cluster, panel (2) shows the PDMF 
obtained between 4 and 9 pc, the PDMF obtained between 9 and 14 pc is shown in panel 
(3), and in panel (4) we show the PDMF corresponding to 14 and 20 pc. 
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Fig. 6. — Variation of the PDMF slope T as a function of the radial distance from the center 
of NGC 346. 
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Fig. 7. — Spectral classification of the stars as derived by iMassey et al.l (119891 ): 03 stars 
are in red, 04 in orange, 05 in yellow, 06 in green, 07 in cyan, 08 in blue and 09 in purple. 
The white open star indicates the center of NGC 346, and the white circle the half mass 
radius, located at 9 pc from the center. 



